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Рассматривается система уравнений газовой динамики горения на основе исходных понятий «нормальная скорость» и «поверх- 
ность горения». Турбулентное пламя описывается одним нелинейным параболическим уравнением. Проведены расчеты взрыв- 
ных процессов в трубах для ряда горючих смесей. Сравнение теоретических данных подлине и времени перехода медленного 
горения в детонацию с экспериментальными результатами показывает их удовлетворительное согласие. 


Изложенные в [1] рассуждения привели к пред- 
ставлению о турбулентном пламени в трубе радиуса 
й как о хаотически блуждающей поверхности с 
фрактальной размерностью сі где ее элементарный 
участок движется относительно газа с нормальной 
скоростью и„ ламинарного пламени. Поверхность 
горения топологически может быть и многосвя- 
зной. Весь газ состоит из двух компонент - свежей 
смеси с массовой долей С (концентрацией) и про- 
дукта горения. Как и С, все термодинамические и 
гидродинамические параметры потока газа по от- 
ношению к гидродинамическому хаосу имеют та- 
кой же смысл средних величин, который они име- 
ют в не турбулентных средах по отношению к хао- 
су молекулярному. Т.е. остающиеся макроскопиче- 
скими временные и пространственные масштабы 
турбулентности много меньше аналогичных мас- 
штабов рассматриваемых ниже процессов. Для 
пламени в таком макроскопическом описании 
можно ввести коэффициент диффузии О [2]: 

0 = 2аи„(1 + В^-У, у = сі, - 2, (1) 

Ип 

содержащей турбулентную пульсацию и ' скорости 
потока газа. Превращение исходной свежей смеси 
в продукт горения описывается выражением для 
скорости химической реакций Ф 

Ф(С)= (1 + 5 и ' уС(\~С)\ (2) 

й й/) 

приходящим на смену закону Аррениуса. Свобод- 
ный параметр В зависит от свойств газовой смеси. 

Система уравнений газовой динамики горения 

Пусть р, р, Т и и - средние по сечению трубы 
давление, плотность, температура и скорость газа. 
В общих случаях без ограничений на скорость дви- 


жения турбулентного пламени, газодинамические 
параметры и концентрация свежей смеси С, могут 
быть найдены решением уравнений [2] 


Ші 

ді 


-чТ= Г — Т-{%*М+Я-ФІС), 
дх у р\ді дх ) с р 

ри + -^~ (ри 2 + р) = ~— ри \и\ 9 
ді дх а 1 1 


др дри 


+ = о, р = рК 8 Т, 2) = 2и„а-К(и', и „ ), (3) 


Ф(С) = —К(и\ и„)-С(1-С) 2 ,К(и', и п )=(\+В^~у 
а и„ 


где 0 - тепловой эффект горения горючей газовой 
смеси с теплоемкостью с р при постоянном давле- 
нии и газовой постоянной К р , у - показатель ади- 
абаты; і,х — время и координата вдоль оси трубы. 

Уравнения (3) не замкнуты: не хватает связи 
между пульсацией скорости и скоростью и потока 
газа. Также система (3) должна быть дополнена 
связью фрактальной размерности с параметром, 
характеризующим турбулентное пламя. Таких па- 
раметров здесь две: число Рейнольдса Не и отноше- 
ние и'/и„. В дальнейшем мы будем пользоваться 
простейшей формулой замыкания [3] и— Ѵя|и|. Что 
касается фрактальной размерности поверхности 
горения, то функциональную зависимость 
4 = 4(Ке,м'/ы„) пока не удается установить. 

На ряде примеров покажем непротиворечи- 
вость ур. (3) ранее установленным теоретическим и 
экспериментальным фактам. Известно [4], что с 
точностью до ошибок эксперимента концентра- 
ционные и тепловые пределы ламинарного горе- 
ния и детонации или совпадают, или мало различа- 
ются. В качестве примера в табл. 1 приведены экс- 
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периментальные данные по концентрационным 
пределам детонации и нормального распростране- 
ния пламени ряда газов в смеси с воздухом [4]. Из 
(3) следует, что эти пределы совпадают. Действи- 
тельно, если и= 0, что имеет место за концентра- 
ционными пределами, то ур. (3) будут описывать 
химически инертную среду, где заведомо не может 
быть ни детонации, ни медленного горения. Таким 
образом, сформулированная модель горения газов 
по концентрационным пределам двух предельных 
режимов - дефлаграции и детонации не противо- 
речит эксперименту (возможность описания дето- 
национного горения показана ниже). 


Таблица 1. Концентрационные пределы детонации и мед- 
ленного горения для ряда горючих газов в смеси с 
воздухом 


Горючее 

Концентрационный предел детонации (и нормаль- 
ного распространения пламени), % 

Нижний 

Верхний 

Этан 

2,87(3,0) 

2,57 (2,20) 

Пропан 

1,98 (1,90) 

3,32 (3,10) 

Бутан 

3,55 (2,25) 

12,2 (12,5) 

Этилен 

7,37(9,50) 

6,18 (8,50) 

Пропилен 

14,7(32,0) 

10,4 (10,8) 


В [5] анализом ур. (3) показано, что в однород- 
ной турбулентной среде (и-соші) пламя распро- 
страняется со скоростью 

и ( = и п (1 + Ви'/и п У, ( 4) 

и что она является собственным значением нели- 
нейной краевой задачи. Последнее обстоятельство 
является одним из важней ш их математических ха- 
рактеристик волн горения [6] . При отсутствии тур- 
булентности скорость пламени в точности равна 
экспериментально наблюдаемому значению: и ( =и я . 
Если интенсивность турбулентности высока, то вы- 
ражение (4) переходит в формулу Вилюнова-Кли- 
мова. То., используемые представления горения, 
основанные на понятиях «поверхность горения» и 
«нормальная скорость пламени» не противоречат 
известным полуэмпирическим и статистическим 
подходам в теории турбулентного горения. Более 
того, если ранее число 5 в формуле Вилюнова-Кли- 
мова присутствовало как формальный параметр, то 
здесь оно имеет явную геометрическую трактовку, 
указывая на меру сложности поверхности горения. 

Расчет перехода медленного горения в детонацию 

Переход медленного горения в детонацию бы- 
строгорящих смесей изучался на основе системы 
(3), в которой первое уравнение переписано в удоб- 
ной для численного решения форме и всем уравне- 
ниям придана дивергентная форма. Переход к без- 
размерным величинам проводился по следующим 
масштабным параметрам 

и, = и П М , р, =р 0 ,р, =р 0 и 2 „, 

К =а/и,„х, = аМ, М = ^р 0 / р 0 и 2 . 


Здесь р 0 , ро - начальные значения давления и 
плотности газовой смеси. Численно решались ура- 
внения в безразмерных переменных: 
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Выражение для коэффициента сопротивления 
Я (формула Никурадзе) взято из [7]. Коэффициент 
кинематической вязкости ѵ считался постоянным. 

Сформулируем начальные условия. При выводе 
выражений для Ф( С) ( 1) и П (2) заведомо предпола- 
гается наличие фронта горения. Это означает, что 
ур. (3) или (5) принципиально не могут описывать 
стадию воспламенения горючей смеси. Начальной 
момент горения должен моделироваться методами 
классической теории горения, изложенной, напри- 
мер, в [6]. Поэтому при решении (5) начальные 
условия необходимо задавать в виде сформировав- 
шегося фронта пламени. 

Потери тепла в стенку трубы в расчетах не учи- 
тываются. Начало взрывного горения происходит в 
центре трубы неограниченной длины. Это эквива- 
лентно зажиганию у закрытой части при адиабати- 
ческих условиях на жестком торце. 

Начальные условия имеют вид: 
р = 1,5, С = 0, и = 0, р = 0,3, для - г < х< +е; (6) 

р = 1, С = 1, и = 0, р = 1 - для других х. 


Константа е подбиралась таким образом, чтобы 
размер области «начального» горения для всех сме- 
сей было равным 3 а. Показатель адиабаты и коэф- 
фициент кинематической вязкости брались равны- 
ми 1,4 и 1,510“ 5 м 2 /с. Радиус трубы о= 1 ,3- 10 2 м. 

Система ур. (5, 6) решалась численно по явной 
2-го порядка точности по пространственной коор- 
динате схеме Годунова-Колгана. 

Моделирование перехода медленного горения в 
детонацию проводилось для смесей 4Н 2 +0 2 , 
2Н 2 +0 2 , Н 2 +0 2 , Н 2 +20,, СН 4 +20 2 , С 2 Н 2 +воздух, 
2С0+0 2 , по которым имеются наиболее полные 
экспериментальные данные, включая условия про- 
ведения опытов. Все начальные данные соответ- 
ствуют нормальным условиям. 

Сравнение расчетной длины К в преддетона- 
ционного участка и времени Г в перехода медленно- 
го горения в детонацию с экспериментальными 
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данными [8-10] показывает их удовлетворительное 
согласие (табл. 2). Основные размерные параметры 
(начальные давление и плотность топливной сме- 
си, тепловой эффект реакции, нормальная ско- 
рость пламени и пр.), определяющие К п и 1 В , отве- 
чают условиям проведения экспериментов. 

Если параметр В подобран для согласования 
расчета с экспериментом, то близость значения по- 
казателя 5 к единице легко обосновать. В экспери- 
ментальных исследованиях факельного турбулент- 
ного горения газовых смесей значение 5«0,7...0,8 
[10]. При сильной турбулентности, имеющей место 
в период ускорения пламени, фрактальная размер- 
ность поверхности горения 5-»1. Но предельно до- 
пустимое значение 5=1 (т.е. ф= 3) реализуется лишь 
при Ке=со. Поэтому здесь и далее 5=0, 8. ..0,9. 


Таблица 2 Сравнение теории с экспериментом по определе- 
нию протяженности преддетонационного участка 
К 0 /2а и времени перехода і 0 


Смесь 

0(5) 

і «о/2 а 

Іо, мс 

теор. 

эксп. 

теор. 

ЭКСП. 

2Н 2 +0 2 

25(0,8) 

45 

38 

1,4 

1,2 

Н 2 +20 2 

15(0,9) 

86 

92 

9,5 

8,3 

4Н 2 +0 2 

30(0,8) 

54 

79 

5,9 

3,1 

Н 2 +0 2 

30 (0,9) 

44 

41 

1,6 

1,8 

СН 4 +20 2 

20(0,9) 

57 

56 

2,1 

2,1 

С 2 Н 2 +воздух 

22(0,9) 

67 

92 

5,7 

6,1 

2С0+0 2 

40(0,9) 

39,6 

37,5 

3,8 

3,0 


Эволюция распределения давления и скорости 
для смеси Н 2 +20, приведена на рис. 1-3, для смеси 
СН 4 +20, - на рис. 4-6. Центру симметрии соот- 
ветствует х=0. На остальных рисунках индексы у 
символов указывают на приведенные в подрису- 
ночных подписях моменты времени. 

При сверхзвуковой скорости движения турбу- 
лентного пламени в его области расположения на- 
блюдаются быстрые пространственные изменения 
всех физических параметров (рис. 1 , б), характери- 
зующих горение. Так как значение приведенной 
скорости химической реакции Ф(С) практически 
целиком определяется турбулентной составляю- 
щей, то положение максимумов на профиле скоро- 
сти и совпадает с максимумами Ф(С). Окрестности 
точек максимумов функции Ф(С) можно считать, 
как и в классической теории ламинарного горения, 
областями локализации турбулентного пламени. 
Они характеризуются высокой температурой и бы- 
стрым изменением концентрации исходного про- 
дукта. На рис. 1, б, видны два следующих друг за 
другом турбулентных пламени. Такая двойствен- 
ность объясняется тем, что за первым фронтом тур- 
булентного пламени формируется волна разреже- 
ния, где скорость газа (как и давление и плотность) 
быстро падает. Это влечет за собой снижение скоро- 
сти химической реакции Ф(С). В результате исход- 
ная смесь не успевает сгореть. Во втором фронте 
пламени снова активизируется скорость горения, 
что вновь сопровождается интенсификацией горе- 
ния, ростом давления, скорости и плотности газа. 




Рис. 1 . Стадия ускорения фронта пламени в смеси Н г +20 2 
для і: а) 7,2, 6) 9,0 мс: 1 ~ р, 2 ~ и, 3 ~ р, 4 ~ С 

Скорость движения второго фронта турбулент- 
ного горения превышает скорость первого. Это об- 
условлено тем, что на процессы горения, согласно 
принятой здесь модели, существенным образом 
сказывается влияние газодинамических факторов. 
Скорости движения фронтов горения зависят от 
характера эволюции (следующих друг за другом) 
ударных волн, которые сами же их генерируют. 

Если скорость пламени меньше скорости звука 
в исходной смеси, то пространственное распреде- 
ление газодинамических параметров и концентра- 
ции С гладкое, без ярко выраженных изломов 
(рис. 1, а). Здесь максимальная скорость газа соста- 
вляет несколько десятков м/с, а скорость турбу- 
лентного пламени относительно газа немногим 
превышает скорость ламинарного горения и„. 

На рис. 1, а, отчетливо видна в точке х&2,6 сла- 
бая волна сжатия впереди фронта пламени, нали- 
чие которой обусловлено начальными данными (6). 

На рис. 2 приведены этапы формирования детона- 
ционного режима горения. Вторая волна турбулент- 
ного горения здесь уже догнала переднюю волну, в ре- 
зультате чего образовалась одна область турбулентно- 
го горения с более высокой скоростью химической 
реакции, большими скачками газодинамических па- 
раметров, чем первые две. После такого слияния волн 
турбулентного горения переход медленного горения в 
детонацию еще больше ускоряется. 
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Рис. 2. Возникновение детонационной волны в смеси 
Н 2 +20 2 для I: 1) 9,4, 2) 9,7, 3) 9,9 мс 


а) 
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Рис. 3. Переход пересжатой детонационной волны в стацио- 
нарную в смеси Н 2 +20 2 


Волнообразный вид распределений давления, 
скорости и плотности за фронтом пламени вызван 
блуждающими возмущениями, возникшими при 
догорании остатков свежей смеси, и частично - при 
слиянии двух первоначальных фронтов турбулент- 
ного горения. Эти волны, взаимодействуя с ударной 
волной и «стенкой» (с центром симметрии), с тече- 
нием времени исчезают. На достаточно больших 
временах, когда завершено формирование стацио- 
нарной детонационной волны, а последняя значи- 
тельно удалена от стенки, блуждающие волны прак- 
тически не наблюдаются. Максимум давления нахо- 
дится всегда немного впереди пламени. Положение 
последней определяется по максимуму скорости хи- 
мической реакции. Образовавшаяся впоследствии 
пересжатая детонационная волна (рис. 3) выходит с 
течением времени (с точностью до ошибок аппрок- 
симации) к режиму, соответствующему на детона- 
ционной адиабате точке Чепмена-Жуге. Экспери- 
ментальное значение скорости детонационной вол- 
ны в смеси Н 2 +20 2 «1927 м/с. Расчетная же скорость 
в используемой модели — 1857 м/с. Также несколько 
заниженной оказалось отношение плотности газа в 
ударной волне р ув к р 0 : 5,6 вместо 6,0. 



0,0 0,5 1,0 1,5 2,0 

х 

Рис. 4. Ускорение слабо турбулентного пламени в смеси 
СНц+20 2 для I: а) 0,79, б) 1, 18мс;1-р,2-и,3-р,4~С 

Качественная картина развития процесса пере- 
хода медленного горения в детонацию для смеси 
СН 4 +20 2 , не отличается от той, что имеет место для 
смеси Н 2 +20 2 . Здесь тоже могут наблюдаться два 
фронта (волны) турбулентного горения (рис. 4, 5). 
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На рис. 4 показано последовательное развитие вто- 
рого фронта, который движется и усиливается бы- 
стрее первого. В нем скорость турбулентного горе- 
ния пока еще относительно низка. Если пик давле- 
ния в первой волне за время от 0,79 (а) до 1,18 (б) 
мс увеличился всего на ед. %, то во второй волне 
давление возросло примерно в 1 ,4 раза от первона- 
чального значения. Эволюционное развитие вто- 
рого фронта турбулентного горения идет столь бы- 
стро, что через еще 0,24 мс (рис. 5, а) амплитуда да- 
вления в нем становится в 1 ,68 раза больше ампли- 
туды давления в первом фронте горения. 



0,0 0,5 1,0 1,5 2,0 


х 



0,0 0,5 1,0 1,5 2,0 


X 

Рис. 5. Фронт турбулентного пламени с двумя очагами горе- 
ния в смеси СН 4 +20 2 для I: а) 1,42, б) 1,71 мс; 1 ~ р, 
2-и,3~р,4-С 

Также сильно различаются и скорости движения 
указанных волн турбулентного горения. Скорость 
первой волны меняется со временем слабо, на рис. 4, 
5 она * 1 100 м/с. Скорость движения второй волны 
быстро растет: если в интервале времени на рис. 4 она 
имеет среднюю скорость «300 м/с, то в промежутке от 
1,18 (рис. 4) до 1,58 мс (рис. 5) составляет «2000 м/с. 

Несмотря на незначительное отличие этой вели- 
чины от скорости стационарной детонационной 
волны 2146 м/с пики давления и плотности в удар- 
ной волне впереди фронта пламени (рис. 5) в 1,5 и 
более раз ниже, чем при детонационном горении. 
Слияние двух волн турбулентного горения в одно 
целое происходит ко времени 1,86 мс. В этот момент 
область максимального давления не имеет остроко- 
нечного пика (рис. 6, о), как при других временах. 
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Рис. 6. Завершающая стадия возникновения детонационной 
волны в смеси СН А +20 г для 1: 1) 1,86, 2) 1,97, 3) 2, 1 мс 

В смеси СН 4 +20, наблюдается переход медлен- 
ного горения к стационарной детонационной вол- 
не через пересжатое состояние. Переход из перес- 
жатой в обычную детонационную волну происхо- 
дит с отрывом от фронта детонации небольшого 
возмущения. Это возмущение, слабо заметное на 
рис. 3, легко обнаруживается на рис. 6, а, как выпу- 
клость на профиле давления позади детонационно- 
го пика. Оно распространяется в противополож- 
ную сторону по отношению к направлению движе- 
ния детонационной волны. Ретонационная волна, 
возникающая в смеси 4Н 2 +0 2 , имеет ярко выра- 
женную форму ударной волны (рис. 7). При прове- 
дении расчетов для этой смеси обнаружено, что ам- 
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плитуда ретонационной волны растет по мере уве- 
личения амплитуды ударной волны, идущей впере- 
ди ускоряющегося фронта турбулентного пламени. 



Рис. 7. Ретонационная волна (х = 0,5) и волна начального 
возмущения (х = 6,4) в смеси 4Н 2 +0 2 . 1 - р, 2 ~ и, 
3 ~ р, 4 ~ С 

Физическая природа появления обратной вол- 
ны заключается в необходимости «уноса» некото- 
рого избытка энергии из пересжатой детонацион- 
ной волны для того, чтобы последняя могла выйти 
к состоянию стационарной детонационной волны. 
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Отрицательное направление движения волны воз- 
мущения, которую называют ретонационной вол- 
ной [8, 10], следует из закона сохранения импульса. 

Выводы 

— Уравнения газовой динамики горения газовых 
смесей допускают формулировку на основе ис- 
ходных понятий «поверхность горения» и «нор- 
мальная скорость» распространения пламени. 

— В такой формулировке уравнений газовой ди- 
намики горения при моделировании турбулент- 
ного пламени отсутствует проблема замыканий 
уравнений типа Фридмана-Келлера, но возни- 
кает другая - проблема нахождения численных 
значений меры Хаусдорфа и фрактальной раз- 
мерности поверхности горения. 

— Разработанная теория дает корректную физиче- 
скую картину перехода медленного горения в 
детонацию, а рассчитанные для газовых смесей 
теоретическое время и протяженность предде- 
тонационного участка удовлетворительно согла- 
суются с экспериментальными данными. Также 
согласие теоретических результатов с экспери- 
ментом имеет место по возникновению перес- 
жатой детонационной и ретонационной волн. 
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